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We investigate dimensional reduction, the property that the critical behavior of a system in the 
presence of quenched disorder in dimension d is the same as that of its pure counterpart in d — 2, and 
its breakdown in the case of the random-field Ising model in which both the interactions and the 
correlations of the disorder are long-ranged, i. e. power-law decaying. To some extent the power-law 
exponents play the role of spatial dimension in a short-range model, which allows us to probe the 
theoretically predicted existence of a nontrivial critical value separating a region where dimensional 
reduction holds from one where it is broken, while still considering the physical dimension d = 3. By 
extending our recently developed approach based on a nonperturbative functional renormalization 
group combined with a supersymmetric formalism, we find that such a critical value indeed exists, 
provided one chooses a specific relation between the decay exponents of the interactions and of the 
disorder correlations. This transition from dimensional reduction to its breakdown should therefore 
be observable in simulations and numerical analyses, if not experimentally. 



PACS numbers: ll.lO.Hi, 75.40.Cx 
I. INTRODUCTION 

The random-field Ising model is a prototypical example 
of a system in which the presence of quenched disorder 
has a drastic effect on the collective behavior. While 
long-range ferromagnetic order is still observed above a 
critical dimension, as in the pure Ising model, the prop- 
erties of the associated critical point are strongly mod- 
ified by the random field. In a renormalization group 
(RG) setting, temperature is irrelevant at the fixed point 
and the long-distance physics is therefore dominated by 
disorder-induced, sample-to-sample fluctuations rather 
than by thermal fluctuations.'^ ^ This leads to anomalous 
scaling relations due to the existence of an additional 
critical exponent associated with temperature, to a shift 
of the lower and upper critical dimensions, a very slow 
spatial decay of the correlation functions at criticality, 
and a very strong slowing down of the dynamics close to 
the critical point that can be described in te rms of an 
unconventional activated dynamic scaling.l^EIIl 

One of the puzzles about the critical behavior of the 
RFIM was the property of dimensional reduction, accord- 
ing to which the behavior in the presence of a random 
field is the same as that of the pure system in two dimen- 
sions less. This property, found to all orders in conven- 
tional perturbation theorj'sMI ^jjd also nonperturbatively 
derived as a consequence of an underlying supersymme- 
try of the model at zero temperatureSffiH^ ■y^,a,s proven to 
be wrong in c? = 3E2I13] pj-om the supersymmetric ap- 
proach it was understood that the failure was related to 
the presence of metastable states, i.e. of multiple ex- 
trema of the bare action (microscopic hamiltonian) in 



the region of interest .'i^l However, no further progress had 
been made. 

We found a resolution of the dimensional-reduction 
puzzle by means of a nonperturbative functional RG (NP- 
FRG) approach, showing that, as in the simpler case 
of a pinned interface in a random environment where a 
perturbative FRG analysis is sufficient ,'^^1'^ breakdown 
of dimensional reduction is related to the appearance 
of a singularity in the functional dependence of the cu- 
mulants of the renormalized disorder, with however the 
singularity becoming too weak to cause a failure of di- 
mensio nal red uction above a nontrivial critical dimension 
d ~ 5.1.1221221 More recently, we also showed that dimen- 
sional reduction breakdown is related to a spontaneous 
breaking of the underlying supersymmetry along the RG 
flow2lH26] g^jjj]^ that it is physically associated with the 
large-scale properties of the avalanches characterizing the 
behavior of the system at zero temperature.^^ 

Whereas the whole description obtained through the 
NP-FRG is consistent and leads to predictions, e.g for 
the critical exponents in d = 3 and d = 4, that are in 
good agreement with computer simulations and ground- 
state numerical studies, directly accessing the properties 
of the RFIM at and around the critical dimension of 5.1 
is not feasible by computer studies (not to mention ex- 
periments!). The goal of the present work is to provide 
a way to get around this problem and to allow for a di- 
rect study of a 3-dimensional system. To this end, we 
consider a RFIM with long-range interactions and long- 
range correlations of the random field. The interest in 
long-range models has a long history and comes from the 
fact that the presence of long-range, power-law decaying. 
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interactions decreases the lower critical dimension of a 
model and that varying the exponent of the power law in 
a fixed dimension allows one to find a spectrum of critical 
behavior that goes from mean-field for truly long-range 
interactions to the absence of transition for short-range 
decay while spanning a continuous range of nonclassical 
behavior in between. In some sense, changing the expo- 
nent of the power law at fixed dimension is like changing 
the dimension in a short-range model. 

To study dimensional reduction and its breakdown, 
one must introduce not only long-range interactions but 
also long-range correlations of the random field. As will 
be explained in more detail below, this is the only way to 
produce a supersymmetry in the field theory at zero tem- 
perature and therefore to possibly generate a dimensional 
reduction property. In addition, the exponents charac- 
terizing the decay of the interactions and of the disorder 
correlations have to be related in a specific manner. The 
problem can then be tackled through an extension of the 
NP-FRG approach combined with the supersymmetric 
formalism, which we have previously developed for the 
short-range RFIM.'^StHl xhe main outcome of the theory 
is that there is a nontrivial critical value of the power 
exponent describing the spatial decay of the interactions 
that separate a domain where dimensional reduction is 
valid (for longer-range interactions) from a domain where 
it is not (for shorter-range interactions). This opens the 
way to a direct check of the transition between the pres- 
ence and the absence of dimensional reduction in the crit- 
ical behavior of the RFIM in the physically and techni- 
cally accessible dimension d — 3. 



II. LONG-RANGE MODEL AND 
SUPERSYMMETRY 

The model that we investigate is the field-theoretical 
version of the RFIM with long-range interactions and 
disorder correlations, with bare action (Hamiltonian) 

S[Lp;h] ^ Ssif] - / h{x)ip{x) = 

J X 

- ^ ^ A(x - vMxMy) -I- ^ I '^'fixf + ^_v{x)^ 
h{x)ip{x) 

(1) 

where = / d!^x and the interaction goes as X{x — y) ~ 
\x — when |a; — 1/| ^ 1, with cr > 0; h{x) is a ran- 

dom source (a random magnetic field) that is taken with a 
Gaussian distribution characterized by a zero mean and 
a variance h{x)h{y) — A(x — y) ~ |a; — y\~'^'^~P'^ when 
la:; — 2/1 >> f . For cr > 2, one obviously recovers a model 
with short-range interactions and a similar reasoning ap- 
plies for p < 0. 

Let us repeat the steps of the Parisi-Sourlas supersym- 
metric construction.^ The critical behavior being con- 



trolled by a zero-temperature fixed pointP^, one can fo- 
cus on the ground-state configuration which is solution 
of the stochastic field equation 



5S[ip] h] 
Sip{x) 



J{x)., 



(2) 



where J is an external source (a magnetic field) conjugate 
to the </? field. When the solution is unique, which is 
precisely the crux of the problem and will be addressed 
later on, the equilibrium (Green's) correlation functions 
of the field are obtained from the generating functional 



h- J 



det 



SipS(p 



X exp / J{x)ip{x) 



(3) 

Because of the assumed uniqueness of the solution, 
the absolute value of the determinant present in the 
right-hand side can be dropped and the functional can 
be built through standard field-theoretical techniques .f^Sl 
One introduces an auxiliary bosonic "response" field (p{x) 
to exponentiate the delta functional and two auxiliary 
fermionic "ghost" fields ip{x) and 'ip{x) to exponentiate 
the determinant. In the resulting form, the average of 
the Gaussian random field can be explicitly performed 
and one obtains 



Z[J,J,K,K]^Zh[J,J,K,K] 

VipViftV^jD^J exp <^ - S'ss [<y3, -0, -0] , -I- 



j (^J{x)ip{x) + ip{x)K{x) + K{x)iIj{x) + J{x)ip{x) 

_ (4) 

where two fermionic sources, K{x),K{x), linearly cou- 
pled to the ghost fields have been introduced and 



Sss = / 'fix) 



Sv{x) Jx Jy 5if{x)6ip{y) 
(f{x)A(x - y)(p{y) . 



(5) 



The (/3-field connected correlation functions of the origi- 
nal problem are obtained from functional derivatives of 
W[J, J, K, K] = log Z[J, J, K, K] with respect to J that 
are further evaluated ior K ^ K = J = 0. 

The next step of the construction is to introduce a 
"superspace" by adding to the d-dimensional Euclidean 
space with coordinates x = {x'^} two anti-commuting 
Grassmann coordinates 9, 9 (satisfying ff^ = (F = 09 + 
99 = O^j^, so that the original field and the auxiliary 
fields can be grouped in a single "superfield" <^{x) — 
ip{x) + 9-il}{x) + ^p{x)9 + 99'f{x), where x = {x,9,9) de- 
notes the coordinates in superspace. A similar procedure 
applies to the sources that can be grouped in a single 
"supersource" J{x). At this stage, we leave unspecified 
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the metric of the superspace. By using the properties of 
the Grassmann variables, Eq. [1] can then be rewritten 
in the following compact form: 



2?$exp -5,3 [$ 



where 



Sssm 



2 4! 
X {-\{x-y)- A(.T - y)dsd^^{y, 



(6) 



(7) 

and where = // dddd and ^ = 

To describe the long-distance physics, one needs infor- 
mation only about the low-momentum behavior of the 
Fourier transform of the long-range functions, A(q) and 
A(g), namely 



A(g) = A(0)-Zlh('z')5 + ' 

A((7) = ALfl.(92)-f +A + -- - , 



(8) 



where the higher-order terms in indicated by the el- 
lipses have been dropped as irrelevant. The above ex- 
pressions have their counterpart in real space in terms of 
fractional derivatives. 

For the short-range model, with Zl_r = Al_r ~ 0, 
terms in the interaction and in the disorder correlation 
can be combined to form a "super-Laplacian" in super- 
space with an appropriately chosen metric.^ To find the 
conditions under which this can be generalized to the 
long-range case, it proves more convenient to work in 
Fourier space for both the Euclidean and the Grassmann 
variables. After introducing 77, f] as the Grassmann "mo- 
menta" conjugate to the coordinates ^, Q (a standard def- 
inition of the Fourier transform in Grassmann space is 
use(P^, we rewrite the action S'ss[<I>] in Eq. Q as 



Mb- 



Zq^~f]Tj (Alr {q^y-^ +A 



(9) 

where = / rf'^g/(27r)'^, 77 = {77,77}, ee JJ dijdrj, and 
= (r/2)$2 + (u/4!)$4 with r = r - A(0). 
Assume now an extension of the short-range case with 
a "supermetric" dx^ = dx^dx^ + CdOdO, with C an 
unknown parameter to be determined, and the associ- 
ated "super-Laplacian" A^s = d^d^ + (4/C)9e9g. The 
squared norm of a "supermomentum" q — {(7,77} is then 
given by = q^ + {A/C)r]fj. It is now straightfor- 
ward to check that the long-range components of the in- 
teraction, ZLRXq^y/"^ , and of the disorder correlation, 
Alr {q^y^"^ , can be combined as a power of q^ if and 
only if 



p = 2 - cr 



(10) 



and if the parameter C of the supermetric is chosen ap- 
propriately. Then, 



ZLR{q^r^ ^ ZLR{qr- {i+m^ii'r' 

-ZLR[qr--rvi'^{e)- 



(11) 



which corresponds to the long-range term in the action 
if one chooses C = ^oZl^I A^^i- 

If one also considers the short-range contributions to 
the interaction and the disorder correlation, an additional 
condition is required, that relates Z and A as 



A 



2A 



LR 



LR. 



Z . 



(12) 



and this can be generalized to include higher powers in 
momenta (see the conclusion). Under the above condi- 
tions, the action S'ss[$] can now be reexpressed as 



C/s($(x)) 



$(-<7,5) Zifl(g2)f + V $(g,^) 

(13) 

which is the generalization to superspace and superfield 
of a if>^ action in Euclidean space with long-range inter- 
action and no disorder. (Note that when Zm and Z 
are different from zero, they can simply be set to 1 by a 
simple rescaling of the fields and momenta.) 

As in the short-range case, the above action is invari- 
ant under a large group of both bosonic and fermionic 
symmetries (the latter symmetries mixing bosonic and 
fermionic fields).!^ Of special importance is the super- 
symmetry associated with the orthosymplectic group 
OSp(2,d)p3 that contains the "superrotations" that pre- 
serve the metric of the superspace. As a result of 
the latter, it can be s hown, both perturbativelj^ and 
nonperturbatively,'i^'iiJ that the superfield theory with 
action 5*83 ['I'] for a Euclidean dimension d reduces to 
the simple field theory with action 5*3 s[(^] in dimension 
d—2. if the superfield theory indeed correctly describes 
the critical behavior of the long-range RFIM this proves 
the dimensional-reduction property. One knows how- 
ever that the Parisi-Sourlas construction breaks down 
when th ere ar e multiple solutions of the stochastic field 
equation.'l^'^ This problem was previously resolved by 
two of us for the short-range RFIlVP^'^ei g^^^j extend 
the proposed formalism to the long-range case below. 

Before presenting the nonperturbative functional RG 
used to describe the long-distance physics, we recall a 
few known predictions of the critical behavior of the long- 
range RFIM. We need to first introduce a few definitions. 
As alluded to in the Introduction, the critical behavior 
of the RFIM, be it short or long range, is controlled by a 
zero-temperature fixed point. The renormalized temper- 
ature is therefore irrelevant and characterized by a criti- 
cal exponent > 0. As a result, the spatial decay of the 
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correlations at criticality are described by two "anoma- 
lous" dimensions instead of one. The "connected" pair 
correlation (Green's) function behaves a^ 



{ip{x)Lp{y)) - {ip{x)){ip{y)) 



T 



\x ~ y\d--^+'l 



(14) 



whereas the "disconnected" one, which survives at zero 
temperature, behaves as 



1 



(15) 



with 2rj> f]> rj and 6 = 2 + rj — fj. 

Consider now the long-range RFIM in d — 3 with p = 
2 — (T, which represents the case of interest in the present 
study. From the results of Ref . [311 , one expects several 
regimes for the critical behavior of the model according 
to the value of the exponent a: 

(i) For a < 1/2, a mean- field regime with classical 
exponents; a ^ 1/2 therefore plays the role of an upper 
critical dimension. 

(ii) For 1/2 < CT < 1, a long-range regime with the 
anomalous dimensions fixed to f] = -q — 2 — a, but non- 
trivial values of the other critcal exponents. (This regime 
corresponds to the long-range exchange and random-field 
correlation regime in Bray's terminology.!^ 

(iii) For a > 1: no phase transition; (7 = 1 therefore 
plays the role of a lower critical dimension. 

Contrary to the generic case studied by Bray, there 
are no other regimes, in particular no short-range regime, 
with the specific conditions p = 2 — a, d = 3. Dimensional 
reduction implies that the critical behavior of the model 
is the same as that of the pure, long-range interaction 
Ising model in d = 1. For this latter system, one also 
expects that the "upper critical" value of a is 1/2 and 
that the "lower critical" one is 1, with no short-range 
regime. Note finally that the long-range RFIM in d = 2 is 
not interesting for the present investigation, as, when the 
underlying superrotational invariance is satisfied, there 
is no range of a for which it displays a nonclassical (non 
mean-field) critical behavior. 



III. NP-FRG FOR THE LONG-RANGE MODEL 

The theoretical approach that we use for investigating 
the critical behavior of the long-range RFIM relies on the 
NP-FRG formalism previously developed for the short- 
range RFIM.iSSlSH It combines four main ingredients: 

(1) A replica or multi-copy formalism, in which the 
permutational symmetry among replicas is explicitly bro- 
ken by introducing linear sources acting independently 
on each copy of the original system. Using expansions 
in the number of unconstrained (or "free" ) replica sums 
then gives us access to the cumulants of the renormalized 
disorder with their full functional dependence, which al- 
lows for the emergence of a nonanalytic behavior in the 
field arguments. ^-^-^^iSSi 



(2) An extension of the parisi-Sourlas superfield con- 
struction in the presence of metastable states. We intro- 
duce a weighting factor exp(— /35), involving an auxiliary 
temperature f3~^, to the solutions of the stochastic field 
equation Eq. ([2| when constructing the generating func- 
tional Z. When f3^^ approaches 0, only the ground state, 
i.e. the state with minimum energy or bare action, con- 
tributes to the functional, as desired. 

(3) An exact functional RG formalism. It is a ver- 
sion of Wilson's continuous RG in which one follows the 
evolution of the "effective average action" , which is the 
generating functional of the 1-particle irreducible (IPI) 
correlation functions. The flow with a running infrared 
(IR) scale k, from the bare action at the microscopic scale 
[k = A) to the full effective action at macroscopic scale 
(/c = 0), is governed by an exact RG equation. 

(4) A nonperturbative super symmetry- compatible ap- 
proximation scheme for the effective average action. It 
involves truncations in the "derivative expansion", i.e. 
the expansion in the number of spatial derivatives of 
the fundamental fields, and in the "expansion in num- 
ber of free replica sums", i.e. the cumulant expansion, 
and leads to a closed set of coupled NP-FRG equations 
that can be solved numerically. 

Steps (1) and (2) above lead, via the usual field- 
theoretical techniques^® and an extension of the deriva- 
tion in section ( [ll| ) , to a superfield theory for an arbitrary 
number n of copies in a curved superspace. The generat- 
ing functional that generalizes Eq. (|6| is then expressed 
as 



a=l ^ a=l-^^ 



(16) 



and the multicopy action is given by 

[Ai«(-a2)-i+^/2 ^Ad^]Mx,0,)-^perm{12)^ , 

where (—9^)" with a a real number is a symbolic no- 
tation describing a fractional Laplacian in Euclidean 
space; its Fourier transform generates a term in (g^)" 
(and for a = 1 one recovers the standard Laplacian). 
In the above equations, we have introduced a super- 
space whose 2-dimensional Grassmannian subspace is 
curved, with the curvature proportional to /?. For in- 
stance, the integral over Grassmannian subspace is de- 
fined as Jg EE / /(I -I- (39e)d9dei^ As discussed in detail 
in Ref. [25], the action in Eq. (17) is still invariant under 



a large group of symmetries and supersymmetries. 

We have applied the NP-FRG formalism to this su- 
perfield theory [step (3) above]. This proceeds by first 
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introducing an infrared (IR) regulator that enforces a 
progressive account of the fluctuations to the bare ac- 
tion, 



2 ^ 



a\X2-, ^) 
'i'a{Xi)Rk{\xi - X2\)<^bix2), 

with two IR cutoff functions Rk and Rk suppressing the 
integration over modes with momentum jgl <C fc (but not 
for those with |g| 3> /cj^^ in the modified fc-dependent 
generating functional 

Zi^\{Ja}]=l ni?$aexp(-5(«[{$4] + 
a=l ^ 

n „ 

E / Ja{x)^a{x)-AS[^\{^a}\ 

„_1 J X 



(19) 



We have chosen the two IR cutoff functions to be related 
through 

Rk{e)^-{'^-^)d,.Mq^). (20) 

The above relation and the form of the regulator ensure 
that all symmetries and supersymmetries of the theory 
are satisfied. This includes the superrotational invariance 
found (only) when the theory is restricted to a single 
copy and to an infinite auxiliary temperature, /3 = 0.^^ 
A specific form for the cutoff function Rk will be given 

below. 

We next introduce the effective average actioUjl^^lMl 
which is obtained from log zjf^ through a (modified) Leg- 
endre transform: 

rr[{$a}] = 

-\ogZi^\{J,}] + Y, [ $„(x)J,(x)-A5f 

(21) 

As already mentioned, it is the generating functional of 
the IPI correlation functions j^^Sl and it is the central quan- 
tity of our NP-FRG approach. Its dependence on the IR 
cutoff k is governed by an exact renormalization-group 
equation (ERGE),I231 



5tri«[{$4] = ^Tr{a,7^.. T'f [{$4]} 



(22) 



where t = \og{k/ K) and the trace involves summing over 
copy indices and integrating over superspace; the modi- 
fied propagator T^^^J^ (xi , ) is the (operator) inverse of 
Qp(^/3)](2) where [ri^^^]'^) is the second func- 

tional derivative of the effective average action with re- 
spect to the superfields ^a{x). 



The effective average action can be expanded in 
increasing number of unrestricted sums over copies, 
which generates an analog of a cumulant expansion for 
the renormalized disorder (more details are found in 
Refs. P2 25 ): 

ri^' [{<^a }] = E ri'^i' [^a] - ' E rl^^ [<f . . <J>^] + • • • (23) 



where (with a pinch of salt, see Refs. [22125] ) . the 
pth cumulant of the renormalized disorder at the scale 
k. Such expansions in increasing number of free sums 
over copies lead to systematic algebraic manipulations 
that allow one to derive a hierarchy of coupled ERGE's 
for the cumulants of the renormalized disorder from the 
ERGE for r0{$a}], Eq. g. 

In Refs. |25)26j . we showedthat the ground state dom- 
inance when /3 — ^ cx) comes with a formal property of the 
random generating functional, which was termed "Grass- 
mannian ultralocality" . This property is also asymptot- 
ically found for finite /3 when fc — ?> (after going to di- 
mensionless quantities, i.e. at the fixed point). When it 
is satisfied, the ERGE's for the cumulants simplify. They 
only involve "ultralocal" parts of the cumulants, e.(?.. 



Je 

ri^2^[<5'l>*2]= / / r,2[$l(^i),<i>2(£2)] 



(24) 



etc, where, in the right-hand sides, T^i, Tk2, only 
depend on the superfields at the explicitly displayed "lo- 
cal" Grassmannian coordinates (on the other hand, the 
dependence on the Euclidean coordinates, which is left 
implicit, still involves derivatives). Ffci, Tk2, are 
then shown to be independent of the auxiliary tempera- 
ture /3~^, and the corresponding ERGE's can be further 
evaluated for physical fields '^aix) — 0a (a^), for su- 
perfields that are uniform in the Grassmann subspace.'^Sl 
For instance, one obtains 



- 7.dt I Pk;X2xA<l}l\{^k2^x2.X:, [01, 0l] " Rk;X2X^) 
^ J X2X3 



(25) 



and 



'ti fe2 [<Pl,(P2\ 



{ Pk;X3X4[4'l]^k3;x'3,.,Xi 



Pk,X3X4 [01, 0l] ^k2l3X4,. [01, 02] 



^Pk,x:>X4 [01,02] 



{^k2Lx, [01,02] - Rk;x,x,) +perm{l2)}, 



(26) 



where perm{12) denotes the expression obtained by per- 
muting 01 and 02, dt is a short-hand notation to indi- 
cate a derivative acting only on the cutoff functions (i.e.. 
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dk = dkRk 5/6Rk + dkRk 5/SRk), and superscripts indi- 
cate functional differentiation with respect to the field 
arguments. The propagators Pk and Pk are defined as 



k2m,(P2 



PkU] = r 



.(2)r 
k.ll 



(27) 



and 



with the higher-order cumulants set to zero. 



(31) 



Pk[(f>i,(t>2] = Pk[(f>i]iTk l\(l)i,(f>2] - Rk)Pk[(t)2] ■ (28) and (p6|, the above ansatz provides 5 coupled flow 



Inserted in the ERGE's for the cumulants, Eqs. (25) 



Finally, we use the fact that the continuous symme- 
tries and supersym metri es of the theory lead to Ward- 
Takahashi identitie^^HSD each running scale k. Im- 
portant ones are those associated with superrotational 
invariance when the multi-copy theory is restricted to a 
one-copy theory by an appropriate choice of the sources 
and the auxiliary temperature has dropped from the 
formalism .12^1 In particular, one then has for uniform field 
configurations, 



2A 



LR 



<jZ 



LR 



(29) 



Note that the relation in Eq. ( 20 ) between the two IR cut 



off functions is precisely of the same form. So long as su- 
perrotational invariance is not spontaneously broken and 
Eq. ( 29 ) remains valid, the ERGE for the first cumulant 
Tki can be closed thanks to the above Ward-Takahashi 
identities [see Eq. (25l]. It exactly reduces to that for 



the effective average action in the pure model in di- 
mension (c[— 2), which entails the dimensional-reduction 
property.E^^H 



IV. SUPERSYMMETRY-COMPATIBLE 
APPROXIMATION SCHEME 

To actually solve the ERGE's describing the critical 
behavior of the long-range RFIM [see Eqs. (25 1 and 



(26)], we use the supersymmetry-compatible nonpertur- 



bative approximation sch eme t hat we have already in- 
troduced in previous work .12^1211 Jt combines a truncation 
in the derivative expansion, which approximate the long- 
distance behavior of the IPI vertices, and a truncation 
in the expansion in cumulants of the renormalized dis- 
order, while ensuring that the Ward-Takahashi identities 
associated with the supersymmetry are not explicitly vio- 
lated. In the present case, the derivative expansion must 
be generalized to account for the long-range spatial de- 
cay of the interactions and the disorder correlations. The 
minimal truncation then reads 



rfei[0] 



-ZkiH^m^mA, 



equations for the 1-copy potential Uk{(t>) that describes 
the thermodynamics of the system, the two field- 
renormalization functions Zk{<f} and ZLfi^k{4')^ the 2- 
copy potential Vfe(0i, (/)2) [from which one obtains the lo- 
cal part of second cumulant of the renormalized random 
field Afe((/)i,(/)2) = 90i502Vfc((/)i,02)] and the strength of 
the long-range component of the second cumulant of the 
random field ls.LR.k{4'i,<l>2)- 

The NP-FRG flow equations are supplemented by an 
initial condition at the microscopic (UV) scale A: = A. It 
is given by the bare action that can be recast as 



(32) 



where, without loss of generality, we have set Z^r — I, 
and 

rA2[0i,02] = ^ |^[0i(a;)(-92)-iW2^2(^) 

+ (/.2(x)(-a2)-i+-/Vi(a;)] + A0i(a;)02(a;) 

with Aiij > 0. To ensure that supersymmetry is not 
violated by the initial action, we choose A — {2/a)Z/S.Lii 
[see Eq. 



The first observation is that due to the long-range na- 
ture of the interaction and disorder correlation, char- 
acterized by a nontrivial exponent cr, the propagators 
Pk[4>] and Pfc [</>!, (/'2] appearing in the flow equations have 
a nonanalytic momentum dependence, even away from 
criticality. By introducing the above ansatz, Eqs. ( 30|3I ), 
in Eqs. (27) and (28), one finds for uniform field config- 
urations. 



{Pkm^q' ^5{q + q')Pk{q^:^ 



(34) 



and 



(Pfc['/'l, 02]),,' = 5{q + q') (^(g2)-l+f [ALRA^l.h) + 

\{cj)id(j)i + 02502)ALj^,fc((/.i, 02)] + Afe(0i, ^^2) - Rk{q^) 



Pk{q^;^i)Pkiq^;(t>2), 



(35) 
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with 



(36) 

The RG flow of Z^R^k and ALR,k is obtained by extract- 
ing from the ERGE of the relevant two-point IPI ver- 
tex, r^^l for the former and r[,^2'' for the latter, that 
part which has the proper singular momentum depen- 
dence (when evaluated for uniform fields). One finds as 
a result that the flow of the two functions Z^R^k and 
^L_R,fc involves field derivatives of themselves in such 
a way that if Z^r^u and l^LR,k are independent of the 
fields at the UV scale, which is indeed the case here [see 
Eqs ( [32p3l )], ZLR.,k and ^LR.k do not flow and remain 
equal to their bare values. This is in line with the con- 
clusions of Ref . [3T] . 

In the present 3-dimensional case, any crossover to 
short-range behavior could only occur for cr > 1: indeed, 
the latter behavior is predicted for a > 2 — rjsR and/or 
a > 2 — 2risR + ?7si?, which, in the present random- 
field system where rjsR « 2r]sR — fjsR « 0.5 in d = 3,l2llll 
only takes place in the region where no phase transi- 
tion is observed anyhow. The whole a range of interest, 
i.e. 1/2 < cr < 1 (see the Introduction), is therefore 
in the long-range regime characterized by fixed anoma- 
lous dimensions given hy rj — fj = 2 — a. The latter 
result can be easily understood by considering the prop- 
agators in Eqs. ( 34|35|36 |. When fc — )• at the critical 
point, Pk^o{(f",4>) red uces to the connected pair corre- 



lation function in Eq. (14 1 while Pk-fo{q , 4>, ((>) reduces 



to the disconnected pair correlation function in Eq. (15 1. 
The low-momentum behavior in Eqs. ( 34|35|36 1 is domi- 
nated by the singular terms, which are both in {q^)^"/"^. 
Comparison with the definitions of the anomalous dimen- 
sions in Eqs. ( 14|15 ) then directly provides the result. 

To cast the NP-FRG flow equations in a dimension- 
less form that allows one to describe the long-distance 
physics associated with the critical point, one must in- 
troduce scaling dimensions that account for the fact that 
the flxed point is at "zero-temperature". Near such a 
fixed point, the renormalized temperature is irrelevant 
and characterized by an exponent > 0, and one has the 
following scaling dimensions (see also section [ll| : 



Uk - fc'^-^ Vk - fc'^-^^ 



(37) 



so that the local component of the second cumulant of 
the renormalized random field scales as k~'^^^~^\ As 
discussed above, 77 = 77 = 2 — cr, so that = 2, in the 
long-range regime considered here (with 1/2 < cr < 1 
in d = 3). Note that, in contrast with the short-range 
RFIM, the equality of -q and fj and a fixed value of the 
temperature exponent 6 — 2 do not necessarily mean 
that dimensional reduction is obeyed. 



Recalling that the long-range functions are not renor- 
malized and using lower-case letters, Uk,Vk,5k,'^, to de- 
note the dimensionless counterparts of t/^, V^, A^,, (/), the 
dimensionless form of the fiow equations can be symbol- 
ically written as 



(38) 



where a prime denotes a derivative with respect to the 
argument. The beta functions themselves depend on w^,, 
Zk, 5k and their derivatives. They are obtained from the 
expressions for the short-range RFIM given in Ref. ^2^] 
after the following replacements: 

Zk{v)y + u'li^) + s{y) ^ + zk{^)y + <(^) + s{y), 
Sk{v>i,'P2) + s'{y) ^ y-^+^^^ + 6k{^i,^2) + s'{y), 

(39) 

where y — /k"^ is the rescaled squared momentum and 
s{y) is the dimensionless form of the IR cutoff function 
that is introduced through 



Rk{q^) - k^s{q^/e) 
Rk{q 



2^ = --ALRd.-Mq^) = --ALRs'iq^lk^) . 
a a 



(40) 

Note that even if one chooses an initial condition where 
the short-range components of the interaction and the 
disorder correlation are zero, Z = A = [which still 



satisfies Eq. (12)], these components are generated along 
the RG flow. 

As the chosen ansatz for the renormalized cumulants 
and for the IR cutoff functions do not explicitly violate 
the Ward-Takahashi identities associated with superro- 
tational invariance, we find the same property as in the 
short-range RFIMpS! so long as the local piece of the cu- 
mulant of the renormalized random field 5k (vi , (^2) is well 
enough behaved when ^ (pi, i.e. does not develop a 
cusp in |(/?2 — </'i|, superrotational invariance is not spon- 
taneously broken along the RG fiow. One then exactly 
finds that Sk{f,y^) — Zk{ip) and that the fiows of Zk{(p) 
and u'i.{tp) are identical to those obtained in the NP-FRG 
treatment of the pure model with long-range interactions 
in dimension d — 2. Two situations may be encountered. 
Either this remains true down to the fc — limit, and the 
fixed point describes a critical behavior with dimensional 
reduction, or a cusp appears at a specific ( "Larkin" ) scale 
fci > along the flow and dimensional reduction must 
be broken. 



V. RESULTS IN THE THREE-DIMENSIONAL 
CASE AND DISCUSSION 

The main goal of this work is to investigate in the 3- 
dimensional long-range RFIM the existence of a critical 
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value <jdr separating a region of parameter a where di- 
mensional reduction holds from a region where it does 
not. We have thus looked for the signature of the ap- 
pearance of a cusp along the RG flow, which corresponds 
to the disappearance of the dimensional-reduction fixed 
point. To this end, we have studied the second derivative 
of (^2) with respect to (pi and Lp2 when evaluated 

in the limit ^ i^i'- 



(41) 



Vl=V2=V 



In the absence of a cusp, Sh 2if) is finite whereas it blows 
up whan a cusp first appears. 

The fiow of 6k.2{^) is simply obtained from the one for 
Sk{<fi,(p2) [see Eq. (38)] by using Eq. (41) and assuming 
that 



Sk{'Pl,^2) = Skfiiip) + {l/2)Sk,2{v)i(pi-ip2) 



(42) 



with ip — (l/2)((y9i + ip2), when (p2 — (fi- The associated 
beta function only depends on Zk, uj^, Sk,o, 5k,2 and their 
(field) derivatives. In this already mentioned, 

<5fc,o(</') = Zk{p>) and the flows of Zk and u'f. are the same 
as in the pure system in dimension d ~ 2 (with d = 3 
here) at the same level of the derivative expansion. 

We have therefore solved the two coupled partial dif- 
ferential equations for Zk and uj, numerically in d = 1, 
for a range of a between 1/2 and 1, and we have taken 
the result as an input to solve the partial differential 
equation for Sk,2i^)- For the reduced cu toff fu nction 
s(y), we have used, as in our previous work ,'221211 ^(y) — 
a(l -I- y/2 + y^/12) exp(— y), where the parameter a has 
been optimized via stability consideration^^Stal] ^^^^^ j-^^^g 
been varied to provide an estimated error bar on our re- 
sults (in practice, stable results are obtained for a wide 
range of a, between 1.5 and 6). 

The outcome of our theoretical investigation is that di- 
mensional reduction between the RFIM with long-range 
interactions and disorder correlations in d = 3 and the 
pure model with long-range interactions in d = 1 is valid 
for a < (Jdr and breaks down for a > ajjj^, with a crit- 
ical value crufi sa 0.72 ± 0.03. Indeed, and as illustrated 
in Fig. [ij 5k,2{p>) stays finite down to fc — below aoR 
whereas it diverges at a finite scale above aoR- We also 
display in Fig. [2] the divergence of the "Larkin" scale at 
which Sk,2{f = 0) — >■ c» when a '^dr- 

We stress that dimensional reduction breaks down 
above <tdr. despite the fact that 77 = 77 = 2 — o- (as in 
the pure model) and 9 = 2. In this case indeed, the cor- 
relation length exponent v differs from its value in the 
pure long-range system in two dimensions less. We have 
not explicitly computed v at this stage, as it requires the 
resolution of the full set of coupled differential equations, 
including a function Sk{pi,(p2) depending on two fields 
and displaying a cusp in \ipi — ip2\, which represents a 
much harder numerical task that we postpone to future 
investigation. However, one can show that the dimen- 
sional reduction fixed point disappears in the presence of 




FIG. 1: NP-FRG flow of SkMv = 0) in the 3-dimensional 
long-range RFIM. The initial conditions at k = \ {i.e., t = 0) 
for u'^.{ip) and Zk{(p) = Skfiif) are taken as equal to the solu- 
tion at the fixed-point. The upper (color online blue) curves 
correspond to a < ann ~ 0.72: a = 0.7,0.71,0.72,0.722 
from top to bottom; Sk,2{0) tends to a finite fixed-point value. 
The lower (color online red) curves correspond to cr > aoR' 
a = 0.8,0.75,0.74,0.73,0.729,0.727,0.725 from left to right; 
they all display a divergence at a finite RG scale \tL\ (which is 
indicated for the three rightmost curves by a vertical dashed 
line). Here, the value of the parameter a of the cutoff function 
is equal to 2. 
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FIG. 2: Divergence of the "Larkin" RG time \tL\ when a — >■ 
CTjtj^. The line is a simple power law in (ctoh — ct)~^ . (note 
that an RG time of 23 is enormous and corresponds to an 
observation length scale 10^*^ times the microscopic one). 



a cusp and that the exponent v in the presence of a cusp 
has additional contributions not present in the dimen- 
sional reduction case: VRPiAiid) 7^ vising{d— 2). On the 
other hand, below a^ji we rigorously find VRFiuid) = 

Vlsing{d — 2). 

According to our recent work, the validity of dimen- 
sional reduction in disordered systems controlled by a 
zero-temperature fixed point is related to the scaling 
properties of the large avalanches observed in the evo- 
lution of the ground state under changes of the exter- 
nal source.'2^ Avalanches are the physical source of the 
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appearance of a nonanalyticity in the cumulants of the 
renormaUzed disorder. Dimensional reduction is bro- 
ken when the avalanche contributions are relevant at the 
fixed point and is valid when avalanches lead to subdom- 
inant effects. In the former case [a > (Jdr), the scal- 
ing dimension d/ of the largest typical system-spanning 
avalanches at criticality is equal to {d + i — t])/2, i.e. 
to (5 -I- (t)/2 in the present study; on the other hand, 
for a < aoR, dj — (5 + (t)/2 — A, where A > is 
the eigenvalue associated with a "cuspy" perturbation 
to the cuspless dimensional-reduction fixed point.^^ The 
eigenvalue A increases as a decreases. It is easily de- 
rived that A = (5 — (t)/2 around the gaussian fixed 
point (for cr < 1/2), so that rfy = 1 for = 1/2 [while 
{5 + a)/2 = 2.75]. 

Finally, we conclude this presentation by stressing the 
relevance to computer studies, either simulations or exact 
ground-state determinations, which represents the main 
motivation of this work. In such studies, one considers a 
lattice, "hard-spin" version of the RFIM, which has the 
following Hamiltonian in the long-range case: 



^ N N 

H — ^ JijSiSj — ^ hiS. 



(43) 



4=1 



goes as \xi — Xj\~'^'^~^''^ at large distance, and the ran- 
dom fields hi are Gaussian distributed with zero mean 
and long-range correlator hihj = A{\xi ~ Xj\) ^ \xi ~ 
|~''^~^+'^\ To ensure that the corresponding generat- 
ing functional satisfies superrotational invariance, the fol- 
lowing relation, which is a variant of the Ward-Takahashi 
equation in Eq. ( 29 ) , should hold: 



— A(r) = -CrJ{r) 
dr 



(44) 



where Si = ±1, the interaction term = J{\xi — Xj\) 



where C parametrizes the relative strength of the dis- 
order compared to the interactions. Then, powerful al- 
gorithms exist to determine the exact ground state of a 
sample in a computer time that scales only as a power of 
the system size."^^ This allows one to investigate large 
systems and to perform efficient finite-size analyses of 
the critical behavior at zero temperature.'2^^ The algo- 
rithms can be extended to describe long-range correlated 
disordeil^ and long-range interactions,!^ and for the cu- 
bic lattice, large system sizes should be attainable. A 
systematic study of the model as a function of the pa- 
rameter <T could therefore provide the first direct inde- 
pendent check of our theoretical predictions concerning 
dimensional reduction and its breakdown. 
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